Heat transfer between a nano-tip and a surface 
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We study quasi-ballistic heat transfer through air between a hot nanometer-scale 
tip and a sample. The hot tip/surface configuration is widely used to perform non- 
intrusive confined heating. Using a Monte-Carlo simulation, we find that the thermal 
conductance reaches 0.8 MW.m _2 K _1 on the surface under the tip and show the 
shape of the heat flux density distribution (nanometer-scale thermal spot). These 
results show that a surface can be efficiently heated locally without contact. The 
temporal resolution of the heat transfer is a few tens of picoseconds. 



Understanding the heat transfer between a hot tip and a substrate is a key 
issue in nanometer-scale devices operating at high temperatures such as scanning thermal 
microscopy (SThM) [l] 2] and also thermally assisted data storage Indeed, several of the 
future methods to write or read information are based on the use of heat transfer between 
an Atomic Force Microscope (AFM) probe and the disk, because a spatial resolution below 
100 nm can be achieved due to the small size of the tip [4|]. Contact heating can be used 
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tip damaging fl. Moreover, in such devices, heat is dissipated over large distances if the 
hot part of the tip is micrometric |6(. A non-contact localized heating through air with a 
moving media could be a solution. Different methods as Joule {e.g. [7]) or Peltier heating 
{e.g. jsl) are possible to heat the extremity of the tip. A quantitative estimation of the heat 
transfer through air clearly appears as a key step to design thermally assisted data storage 
techniques. In this paper we report a calculation of the heat transfer between the apex of a 
tip and a surface. We discuss the leading mechanism, the spatial and temporal resolution. 



ting the substrate [3|. But the tip-sample contact is then responsible of 
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rot tip and a cold sample are not in contact, heat can be exchanged through 
8] or radiation [9( . The radiation contribution can be evaluated using the model 
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FIG. 1: A skip of the tip. Typical lengths used are 1=40 nm and h=20 am. The distance d is in 
the nanometric range. 

I — I 

described in ref.[9|. The heat conduction cannot be modelled by classical Fourier diffusion 
because the air mean free path (MFP) is on the order of the average distance between the 
tip and the sample. Pure ballistic transfer is neither valid : this problem has to be solved 
in the intermediate regime. 3D effects have also to be included to observe the impact of the 
tip geometry. We therefore use a Monte-Carlo approach associated with the linear response 
theory in order to directly solve space and time dependent problems. This model also allows 
to study the time dynamics of heat conduction through the air gap. 

We consider the extremity of a pyramidal AFM tip with a square base of I x / with 
Z=40 nm and a height h—20 nm. The distance d between the tip and the surface is in the 
nanometric range. We give the thermal characteristics of heat diffusion in the tip. In the 
stationary regime, heat confinement remains realistic since heat diffusion in the tip generates 
a significant temperature decay over a few tens of nanometers. In the transient regime, the 
heat diffusion in the tip has a characteristic time (h 2 /a, where a is the diffusivity of the 
material) on the order of 1 ns. 
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We now describe how the radiation is evaluated. The radiation between a sphere with 
same volume as the heated zone of the tip and the surface was calculated for different 
materials. For a distance 20 nm between the center of the sphere at 800 K and the surface 
of same material at 300 K , we found a flux of 5 10~ 17 W for Si and 1 10~ 14 W for silver. 

n 

As discussed in ref. [9[, this flux can be increased if both the tip and the surface have an 
infrared active optical phonon resonance. For glass, we found 9 10 -9 W. 

We now turn to the description of the conduction model. Its basis consists in calculating 
the time evolution of the local heat flux on the sample < (r, t, T tip ) > due to an impulse 
of molecular flux leaving the tip. The heat flux response q t h (r, t) to a given molecular flux 
excitation q^f (r) is then given by 

qth (r, *)=/'< q& (r, t - r, T Up (r)) > S tip .q% (r) dr (1) 

J 

where < > plays the role of a susceptibility. S tip is the area of the heated part of the tip. 

Let us now discuss how the susceptibility is computed using a Monte-Carlo approach. The 
heat transfer is due to air molecules which are in first approximation nitrogen molecules. A 
molecule is considered to carry the translational and rotational energy I&bT where fee is 
the Boltzmann's constant (reference 111]). A cold molecule is heated when impacting the 
hot tip and then flies with a velocity chosen according to a Boltzmann's law at T tip . The 
emission laws for velocity v and direction described by the cylindrical angles 9 and ip are 
derived from the equilibrium molecular flux leaving the probe at T tip 

xfl Q cos9sm6d6f*Z d<f. (2) 

where m is the diatomic molecular mass. nu p is the number of molecules per unit volume 
derived from the condition of null molecular flux at the tip surface. The sample and gas 
are at the ambiant temperature T a . Neglecting the increase of T a due to the presence of 
the tip, the incident molecular flux is given by Eq. (2) where n t i P is replaced by the number 
of molecules per unit volume in the gas denoted n g and T tip is replaced by T a . A molecule 
leaving the tip flies until undergoing a collision with another molecule or with the sample 
surface. The path length between two collisions is computed according to an exponential 
law with a characteristic decav length given bv the MFP A = * which is 55 nm at 
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atmospheric pressure. R pot = 0.2075 nm [lOj denotes the molecular radius. The velocity of 



4 



the collision partners are computed with Boltzmann's law taken at ambient temperature T a . 
It allows to define a new velocity and a new direction after the collision for the test-molecule 
according to the Very Hard Sphere model [h]]. A molecule may undergo several collisions 
before it reaches the sample surface. Hot molecules which remain in the gas are discarded 
and the contribution of collision partners to the tip-sample heat flux is neglected. Finally, 
the thermal flux received by an element Sx x 8y of the surface is 

where N is the number of emitted molecules and E (r, t) is the net energy transfered by 
the molecules emitted by the tip and falling in the area 5x5y during the time interval St. 
Each molecule arrives with an energy determined by the Monte-Carlo simulation and leaves 
the area with the mean energy T a . We assume therefore that the kinetic energy is fully 
absorbed Those molecules are considered to leave the tip at the same initial time and 
the arrival time t is given by their flight time. The calculated susceptibility is a statistical 
average over iV~4x 10 8 molecules. At T tip = 800 K, the mean power given by the tip to 
air molecules is about 90 nW. 

Figures 2(a) and 2(b) represent the heat flux density deposited on the sample surface 
when the tip is in contact (top) and when d = 10 nm (bottom) in a steady-state analysis. It 
is obtained by temporal integration of the susceptibility. The characteristic size (FWHM) 
of the heated zone in contact is a square of edge 35 nm. This corresponds to the best spatial 
resolutions for SThM and also to disk densities up to Tbit.in -2 . The geometry of the heat 
flux density reproduces the square shape of the tip projection in the case of contact. For 
heights larger than d = 10 nm, the imprint of the tip is lost and the deposited energy becomes 
axisymmetric. This shows that the shape of the tip no longer influences the heat transfer. 
This behavior is still observed if we do not include collisions in the model. It follows that it is 
due to the isotropy of the velocity distribution of emitted molecules. Previous models based 
on vertical projection might not reproduce correctly this phenomenon under d = 10 nm 
Levels of heat flux densities obtained for tip temperature of 800 K are higher than 45 MWm -2 
in contact. The values that we found correspond to a total power on the order of 10 nW 
much larger than the radiation contribution when using silicon. The temperature of a 10 nm 
thick film of insulator material (thermal conductivity k=l Wm _1 K _1 ) would be increased 
by about 80 K in 1 fis in contact and by 25 K when d — 10 nm. This indicates that 
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FIG. 2: Heat flux density on the substrate, (a) d = nm (b) d = 10 nm (c) Heat flux density along 
the section shown in (a)(b). Dashed lines show results when collisions are not taken into account. 
The inset shows the normalized heat flux density. 
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FIG. 3: Temporal evolution of the susceptibility : power given to a s = 15 nm square of the 
substrate for a I = 40 nm and h = 20 nm tip. Dashed lines represent pure ballistic calculations. 
The inset shows the full power on the sample normalized by the pure ballistic flux when the tip is 
retracted. This illustrates the decay of the flux due to collisions. 

heating through air with a moving media is feasible. On Figure 2(c), we report the radial 
distribution of the heat flux density. In contact, a linear decrease is found until 35 nm. As 
seen in the inset, the FWHM of the thermal spot increases by more than 10 nm when d 
increases from d = to 20 nm. The comparison of dashed lines and plain lines shows that 
the effect of the collisions is to further reduce the heat flux. 

Figure 3 represents the evolution of the susceptibility [3| < qf h (r, t, T) > integrated 
over a square surface of side s — 15 nm under the tip with I = 40 nm and h = 20 nm 
at various heights d. The minimum heating time is equal to the time width of the flux 
response. This is on the order of 10 ps. This is not a limiting factor for writing applications 
since this time is much smaller than the conduction time scale in the tip. As shown in the 
inset, the power given to the sample surface at d — 30 nm is half of what is found in a pure 
ballistic calculation and 65% for a square of edge 25 nm. Although the susceptibilities are 
not qualitatively different, it is necessary to take into account effects of collisions. Note that 
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the collisions are responsibles for the loss of 70% of the power in comparison to the ballistic 
one when the height d is equal to the mean free path. 

In summary, we have studied the heat transfer through air between the hot extremity of a 
tip and a surface. The hot tip/surface configuration is widely used to heat locally a surface. 
We find that for silicon, air conduction dominates the heat transfer. The values of the heat 
flux density reach 0.5 MWm _2 K _1 at a distance of 35 nm from the axis of the tip. The time 
scale of the conduction heat transfer is on the order of tens of picoseconds. These figures 
show that heat transfer through the air at distances on the order of 10 nm is sufficient to 
allow non-contact scanning thermal microscopy or thermally assisted data storage. 
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